New definition of the discrete 
fractional Fourier transform 



cn 

Nicolae Cotfas and Daniela Dragoman 

O . 

£NJ , University of Bucharest, Physics Department, 

£ ■ P.O. Box MG-11, 077125 Bucharest, Romania 

E-mail: ncotfas@yahoo.com, danieladragoman@yahoo.com 



c3 



Abstract. The definition of the continuous fractional Fourier transform is based on 
the use of a system of eigenfunctions of the continuous Fourier transform, namely, 
i-Ch ' the Hermite-Gaussian functions. A discrete fractional Fourier transform completely 

analogous to the continuous fractional Fourier transform can be defined by using an 
adequate discrete counterpart of the Hermite-Gaussian functions. The finite frame 
C$ ' quantization is a discrete counterpart of the coherent state quantization. It is known 

that the quantum harmonic oscillator can be defined in terms of the coherent state 
quantization. By using the finite frame quantization, we obtain a discrete counterpart 
of the quantum harmonic oscillator. The eigenfunctions of the finite oscillator obtained 
^ , in this way represent a discrete counterpart of the Hermite-Gaussian functions adequate 

for a new definition of the discrete fractional Fourier transform. This new definition is 
compared with the main existing definitions of the discrete fractional Fourier transform. 

q 
o 

1 1. Introduction 

The fractional Fourier transform (FrFT), denoted by J 7 ", represents a generalization of 
the Fourier transform J 7 directly related to the evolution operator of the harmonic 
oscillator [9]. The discrete counterpart F a , called the discrete fractional Fourier 
transform (DFrFT), is a generalization of the usual finite Fourier transform F, and 
has important applications in optics and signal processing [191 EHj • 

It is well-known that the eigenvalues of the Fourier transform J 7 are 1, — i, —1, i and 
the corresponding eigenspaces are infinite-dimensional. If the real number a is not an 
integer then the numbers l a , (— i) a , (—1)°, i Q are not uniquely determined. Therefore, 
there exists an infinite number of linear operators which can be regarded as the a-th 
power of J 7 . Fortunately, among these operators there exists a privileged one, and this 
allows a natural choice for J 7 ". It is defined by using the Hermite-Gaussian functions 
as eigenfunctions of J 7 " and a particular choice for the corresponding eigenvalues. 

In the case of the finite Fourier transform F we have a similar problem, and a 
privileged choice for F Q is obtained by following the analogy with the continuous case 
[20| |2T] . The best known choice is based on the use of a discrete counterpart of the 
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harmonic oscillator defined in terms of finite-difference operators and leading to the 
Harper functions [H |3] as a discrete counterpart of the Hermite- Gaussian functions. 

We show that a better discrete counterpart of the harmonic oscillator can be 
obtained by using a discrete counterpart of the coherent state quantization (2J [TTJ [16] , 
which we call the finite frame quantization [Hj. The corresponding eigenfunctions are 
at the same time eigenfunctions of F and approximate the Hermite-Gaussian functions 
better than the Harper functions. It is well-known that better Hermite- Gaussian-like 
discrete functions can be obtained numerically by starting from initial given ones and 
using refinement techniques (solving series of constrained minimization problems [H] 
or a modified Gramm-Schmidt computer orthonormalization [201 I2~T]). We are mainly 
interested in versions of DFrFT useful in quantum mechanics, based on functions with 
physical meaning. Therefore, it is desirable to avoid to use the refinement techniques. 

2. Fractional Fourier transform 



The Fourier transform of a function ifj:M. — >C belonging to fl L 2 (R) is 

1 



C. 



i/j(x') dx' . 



'2ix j^oq 

This transformation can be extended to the Fourier- Plancherel transform on L 2 
is a unitary transformation, 

1 



J 7-1 



.T ) 



e lxx il)(x') dx' 



'2n j _oo 
and satisfies the relations 

J r2 ['i/'](x) = ip(—x) , J-" 4 [^](x) = ip(x) , for any xE) 

The normalized Hermite-Gaussian functions 

\l/ m : R — > R, V m (x) = 1 ^ H m (x) e~^ 2 

y/m\ 2 m v /7T 

are eigenfunctions of J 7 

^ m ] = (-i) m * m , for any me {0, 1, 2, ...} 
and {ty m }m=o i s a complete orthonormal set in L 2 (R). We have 



m=0 



(X)= E(*m^)^[*m](s) 



m=0 



£ H) m *m(x) JZo ^m(x') V(X') dx' 
m=0 



oo 
oo 



E e-f m V m (x) * m (x') 

m=0 



ip(x') dx'. 



For aGl, the ath-order continuous fractional Fourier transform is defined as 



E 

771=0 



e~ ma ^ m (x) m m (x') 



ip(x') dx'. 



(1) 
. It 

(2) 



(3) 



(4) 



(5) 



K a (x, x') = < 



Hm m = ( m + J ) (9) 
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It is periodic, J 7Q + 4 = J 70 , and the definition can be written as 

/oo 
K a (x,x')<ip(x')dx' (6) 
-00 

where 

8(x — x') for a = 

^■(a-signa) / (x 2 +X /2 ) COS ^ 0-2xx' \ „ „ , , _ /-n. 

; exp i „ . J 1 tor 0< a <2 (7) 

^vrlsinfal ^\ 2 sin fa y II W 

<5(x + x') for a = 2. 

The functions \I/ m are at the same time eigenf unctions of the Hamiltonian 

H = -- + -x 2 (8) 
2 dx 2 2 v ; 

of the harmonic oscillator corresponding to distinct eigenvalues 

1 
2 

From the relation 

^ a [* m ]=e-^ mQ ^ m = e-^ Q( ^^^ m for any m£ {0, 1, 2, ...} 
it follows the equality [H] 

JT« = e -f a(H-l) 

whence 

e _ it H = e -i| ^ (1Q) 

The last relation shows that FrFT can be viewed as the evolution operator in the case 
of the harmonic oscillator. 

3. Discrete fractional Fourier transform 

In this paper we consider only the case of a Hilbert space % of odd dimension d — 2s+l. 
The space "H, isomorphic to the standard Hilbert space C d , can be identified with the 
space of the functions 

if> : Z d — > C 

that is, with the space of all the <i-periodic functions 

(f : Z — >C with tp(n+d)=(p(n) for any n6Z. 

A o?-periodic function ip : Z — > C is well determined by its restriction to {0, 1, d— 1} 
or {— s, —s + 1, s — 1, s}, and the scalar product of two functions can be written as 

d-l s 

(V,^) = ^2<p(n)iJ;(n) or (ip,ip) = ^ tp(n)>ip(ri). (11) 

n=0 n=—s 

The set {£j}? = _ a , where Sj : Z — )• C is the ci-periodic function satisfying the relation 
Sj{n) = 5j n for ne{- s, — s+1, s — 1, s} 
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is the canonical ortho normal basis of %. Since 

the definitions of the finite Fourier transform 

F : U — >U and of its adjoint F + : H — > H 
can be written as 

FbP) = 4j E e"^ fc " <p(n), F+[ip] (k) = -L E W (12) 

n=0 n=0 

or 

FMW = ^ E e-T^(n), F+[^](Ar) = 4 I E e^^(n). (13) 

The finite Fourier transform is a unitary transformation, that is, 
pp+ _ p+p = j 

where I : % — > H, Ix = x is the identity operator. For any (peH we have 

F 2 [(p](k) = ip(—k), for any keZ (14) 

and therefore, 

F 4 = I and F" 1 = F + = F 3 . (15) 

For d>4 the eigenvalues of F are 1, — i, —1, i, that is, (— i)°, (— i) 1 , (— i) 2 , (— i) 3 . 
The operators Po, Pi, P2, P3 : H — > "H, where 

Pm = \ E (-i)" fcm F fc = \ E ef tm F fc (16) 

k=0 k=0 

that is, 

P = ±(I + F + F 2 + F 3 ) P 2 = ±(I-F + F 2 -F 3 ) 

Pi = |(I+iF-F 2 - iF 3 ) P 3 = f (I-iF - F 2 +iF 3 ) 

are projectors 

P+=P m , P 2 m = P m for any me {0,1,2,3} 

orthogonal 

P m P^ = for m ^ £ 
and satisfy the resolution of identity 

I = P + Pi + P 2 + P 3 . 

They are the orthogonal projectors corresponding to the eigenspaces of F 

FP m = (-i) m P m for any me {0,1, 2, 3} (18) 
and F admits the spectral representation 

3 3 

F = P - iPi - P 2 + iP 3 = Hr Pm = e ' fm P - ( 19 ) 

m=0 m=0 



(17) 



New definition of the discrete fractional Fourier transform 5 

Any non-null vector from % m = P m ("H) is an eigenvector corresponding to (— i) m . 
Consider for each m G {0,1,2,3} an orthonormal basis {r) mi i, r) m> 2,..., Vm,d m } in W m . 
The system of vectors 

{r)o,i,r)o,2, m,d , 1)1,2, ■ i)i,di, V2,l, 1)2,2, r)2,d2, V3,U V3,2, ■■■,V3,d 3 } 
is an orthonormal basis in %, and the matrix of F in this basis is a diagonal matrix 
F = diag(l, 1, -i, -i , -1, -1 , i, ...,i) 

do times d\ times di times d$ times 

If a G R is not an integer, the a-th power of a complex number z, namely, 

Z a = e Q ( ln l ;2 l+ i ( ar g( 2 )+2fe7r)) ^20) 

depending on a parameter fceZ, is not uniquely determined. For example, 

i a a 2kirai ( \\a _(— ^+2kn)ai 

(—\\ a — e (7T+2A;7r)ai ja _ e (f +2fc7r)ai V ix ; 

Each of the diagonal matrices of the form 

f q = dia g ( r, r, (-i) Q , (-i) a , i Q , n 

represents [I] the matrix in the considered basis of an a-th power of F. In the above 
expression l a ,...,l a may be different, (— i) Q , (— i) a may also be different, etc. The 
particular choice 

F£ = diag( e-T»..., e -T», e~" a , e~"°, e~^ a , ...,e~^ a ) 



I times di times efo times 1^3 times 



corresponds to the operator 

3 33 



Ff = E e"f mQ P m = i £ E e? (4 - a)m F 



m=0 " fc=0m=0 



(22) 

V* p 2 ?^-") sin7r ( fc - Q ) F fc 

e 4 4 sin^(fc-a) r 

fc=0 4 ^ ; 



called the standard discrete fractional Fourier transform 

Among the elements of H there is the real d-periodic function [T7J [18] [23] 

oo 

Q K :Z-^R, g H ( n )= £ e ^^ 2 (23) 

l=-oo 

depending on a parameter k G (0, oo), we call a /mite Gaussian [6]. It is obtained by 
using a Zak type transform [32] , and is a discrete version of the continuous Gaussian 

f K :R-^R, f K (x) = e"f* 2 (24) 

because 

Q K {n) ~ /«(n-v/27r/d), for any n G {-s, -s + 1, s-1, s}. 
For example (see Figured]), in the case d = 21, 

f 0.00239624 for « = 1/3 
1.37592 • 10~ 8 for « = 1 (25) 
3.05311 • 10~ 16 for « = 3. 
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Figure 1. Gaussians A/3, fx, A and the corresponding finite Gaussians 01/3, 0i, 03 
in the case d = 21 (the unit of measure used in the discrete case is ^J2n/d). 



The function g K can be expressed in terms of the Jacobi function 63 as 

= 7S e *(l>&) (26) 
and one can prove [6j [23] the relation 

F[fl K ] = ^ i (27) 

representing a finite countepart of 

f*W = 7*h®- (28) 
The real part as well the imaginary part of F^[g re ] are linear combinations of g K and Qi 

(29) 

The standard discrete fractional Fourier transform is a simple choice, but it does 
not approximate the FrFT in the limit d— >oo (see, for example, the figures from [2T]). 
This property being the major motivation for defining DFrFT, the standard discrete 
fractional Fourier transform is less used in applications. However, based on Weyl's 
formulation of quantum mechanics, generalized by Schwinger, the standard discrete 
fractional Fourier transform could be useful in a self-consistent theoretical framework 



F£[ 0K ] = \ (1 + e— ) Q K + 5^. (1 - e— ; ) 
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describing quantum physical systems characterised by a finite-dimensional space, when 
the continuous limit is not an issue. 

On the other hand, a discrete fractional Fourier transform completely analogous to 
the continuous FrFT can be defined by using an adequate discrete counterpart of the 
Hermite-Gaussian functions and a similar choice for the a-th power of the corresponding 
eigenvalues. In the next subsections we present three alternative known definitions. 

3.1. Definition of an approximate DFrFT based on sampled harmonic oscillator (Fq) 

The most straightforward way to obtain a discrete counterpart of the Hermite-Gaussian 
functions is the sampling method. The functions { r d m }ti=o^ where "& m : Z — y R is the 
c?-periodic function satisfying the relation 

= C m * m for any nE {-s, -s+1, s-1, s} (30) 

with C m a normalizing constant, are neither quite orthogonal not quite eigenfunctions 
of F. Nevertheless, the difference is small [30], and the transformation 

where 



FSM(*)= E 



d-1 



V>(n) (31) 



_m=Q 

may be used as an approximate DFrFT. The operator Fq is not a true a-th power of 
F. An improved version of this definition was obtained by Pei et al. [20l [2Tj . 

3.2. Definition of an approximate DFrFT based on Mehta eigenbasis (F^) 

The Mehta's functions {// m }^ , where pi m : Z — y R is the c?-periodic function [18] 

Hm(n) = c m E V m [(n + id) Jf (32) 

i=-oo V V / 

with c m a normalizing constant, are eigenfunctions of the discrete Fourier transform 

F[ yUm ] = (-i) m /i m . (33) 

The functions fi m are not quite orthogonal, but the transformation 
U^U:i>^F a M [i;] 

where 



F a M [m) = E 



d-l 



tp(n) (34) 



J^e ^ n m {k) n m {n) 

_m=0 

may be used as an approximate DFrFT. Again, the operator F^ f is not a true a-th 
power of F, but the differnce is rather small |30j . 
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3.3. Definition of an exact DFrFT based on Harper functions (Fg) 

In this subsection we present the definition proposed by Candan, Kutay and Ozaktas 
[3] and based on some results obtained by Pei and Yeh [20] , Dickinson and Steiglitz [10] . 
By denoting T> = j-, we get 

whence 

TH = FLT. (36) 

This explains why the eigenfunctions \l/ m of H, corresponding to distinct eigenvalues, 
are at the same time eigenfunctions of J 7 . The finite-difference operator I) 2 , where 

&<p(x) = ^ + £) ~ M * ] + V{X ~ £) (37) 

e 2 

is an approximation of V 2 , and 

^,^9^ 1 2(cos(27rex) — 1) . . 

e 2 

The finite-difference Hamiltonian 

H = V 2 + FV^- 1 (39) 
commutes with the finite Fourier transform 

FH = HF. (40) 

For x = ne and e = 1/ y/d the equation H</? = \(p becomes 

cp[n+l] - 2<p[n] + <f[n-l] + 2 (cos ^ - l) cp[n] = \<p[n\ (41) 

where ip[n] := (p(n/\/d). Since the eigenvalues of H are distinct for odd d (see Figure 
|2]), the eigenfunctions of H are at the same time eigenfunctions of F. The rf-periodic 
normalized eigenfunctions h m of H, considered in the increasing order of the number of 
sign alternations, are called Harper functions [TJ [30] . The transformation 

where 



mm) = E 



ip{n) (42) 



- d-i 

J2^ 4ma hm(k) h m {n) 

_m=0 

is called the discrete fractional Fourier transform. The operator F^- is an a-th power 
of F. 



New definition of the discrete fractional Fourier transform 



9 



4. Frames and the frame quantization 

In this section we present (by following [TT]) a generalization of the Berezin-Klauder- 
Toeplitz quantization [21 [16] , the notion of quantization being considered here in a wide 
sense. It contains as particular cases the coherent state quantization and the finite frame 
quantization. To simplify, a quantization can be regarded as a procedure that associates 
to an algebra of classical observables an algebra of quantum observables. 

We start from an observation set X of data or parameters equipped with a measure 
\i and consider the Hilbert space of real or complex square-integrable functions 



L 2 (X,/x) = <^ <p:X — ► C 



i x 



\Lp(x)\ 2 n(dx) <oo . (43) 



If X is a finite set with the discrete measure n({x}) = 1, the previous relation becomes 



L 2 (X 7/ u) = { ip:X — ► C 



Ei^)i 2<o ° ( 44 ) 



x&X 

and, evidently, L 2 (X, /i) is the space of all the functions ip : X — > C. 

Then we choose an orthonormal system {4>j}j^j C L 2 (X,fi) such that 

Af(x) = ^2\4>j(x)\ 2 < oo for any xeX (45) 

and a Hilbert space K. admitting an orthonormal basis of the form {|ej)} je j. The set 
{\%)}xex, where 

VW(x) ffj 

satisfies the relation 

/ fi(dx) M(x) \x)(x\ = V" \ej)(e k \ / <f>j(x) <j) k (x) fi(dx) = I K 

that is, the resolution of identity 

fi(dx) M{x) \x) (x\ = I K . (47) 



x 



It is called a frame for K,. Generally, the vectors \x) are not linearly independent. 

A classical observable is a function / : X — > C having specific properties with 
respect to some supplementary structure associated to X, like topology, geometry, 
etc. The quantization with respect to the frame {|x)} :re x is defined as the linear map 
f \-> Af, where Af is the linear operator 



Af.K^K, A f = fi(dx) f{x)Af{x) \x)(x\. (48) 

Jx 

Particularly, in the case f(x) = l we have Af = I/c. The function 

Af.X — >C, A f (x) = (x\A f \x) (49) 

whose values are the mean values of A f in the states \x), is called the covariant or lower 
symbol of Af. The function / is named contravariant or upper symbol of Af. 
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4-1- Quantum harmonic oscillator 

Consider the set X = C, identified with the phase space for the motion of a particle on 
the line by using the one-to-one mapping 

R 2 ^C: (q,p)^^=(q + ip) (50) 

and equipped with the measure 

li(dz) = — d 2 z, where d 2 z = dfRzz d3mz . 

71 

The set of functions {(f) n } n eN, where 

<P n : C — ► C, <j) n {z) = e"^ 4= ( 51 ) 



is an orthonormal system. Consider /C = L 2 (IR) and the orthonormal basis {|n)} n6 N, 
where \n) is the Hermite-Gauss function ^ n in Dirac notation. We have 

oo 

(n\k) = 5 n k and \ n )( n \ — I- (52) 

n=0 

The elements of the corresponding frame {|z)} z6 c, where 



are the Schrodinger-Klauder-Glauber-Sudarshan states , or simply the standard coherent 
states [TTJ E2]. They satisfy, among many others, the relations 

(z\z) = 1, - f \z)(z\d 2 z = L (54) 



7T 

The quantization with respect to the system of coherent states {|z)}>ec> called coherent 
state quantization, is defined as the linear map / H- Af, where 

/ : C — ► C 

is a classical observable defined on the phase space and 

A f = - f f(z)\z)(z\d 2 z (55) 

the corresponding linear operator. In the case of the function f(z) = 1 we get A\ = I. 
The linear operator corresponding to the function f(z) = z is 

A z = - / z\z)(z\d 2 z = V Vn + 1 |n)(n + l| (56) 

JV n=0 

that is, the lowering operator 

a\n) = \fn\n — 1) (57) 
having the coherent states \z) as eigenstates 

a\z) = z\z) for any zfC. (58) 
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The linear operator corresponding to f{z) = z is 



A- z = - 

7T 



- / z\z)(z\d 2 z = J^Vn+l \n+l)(n\ (59) 

K ^ C n=0 



that is, the raising operator 

a+\ n ) = Vn + 1 \n + l). (60) 
In view of the linearity of quantization 

A ^{z+z) = 75(« + « + )=9 



(61) 



iV2 

are the usual coordinate and momentum operators. We have \z\ 2 = \{p 2 + q 2 ), but 

00 11 
M = + l n >H = 2^ + ^ + 2 l (62) 

n=0 

that is, A\ z \2 is the Hamiltonian of a translated harmonic oscillator. The eigenvalues of 
A\ z \2 are 1, 2, 3, ... (see Figure[2]) and the corresponding eigenfunctions are the Hermite- 
Gaussian functions \l/ n . 

4-2. Quantum finite- dimensional oscillator 

Consider the finite set X = 7L d x Z^ equipped with the counting measure 

/i({(a, /3)}) = 1, for all (a, /3) G Z d x Z d 
and the normalized finite Gaussian 

g : Z d — > R, g(n) = ; =. (63) 



EL-,(0i(^)) 2 

In the Hilbert space of all the functions tp : Z^ x Z^ — )• C, the set {0j}j = _ s , where 

^. :Z d xZ d ^C, p) = -j= e"^ ft ' g(a+j) (64) 

is an orthonormal system. Consider K, = % and the canonical basis {|£j)}i=_ s satisfying 

s 

(sM = Sjt and |ej><£j| = I. (65) 

j=-s 

The system of vectors {\a, « = _ s , where 

|a,/3)= $>^'g(a+j)k;> (66) 



J=-S 

is a finite frame EJ H3] in H 



1 s 

(a,/3|a,/3) = 1 and - ^ |a, /3| = I. (67) 



d 
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H 



Figure 2. Eigenvalues of Q, A q , H. Ai^^j (in the case d = 21) and of A\ z \2 . 



The finite frame {\a, (3)} s a g=- s is Fourier invariant (up to some phase factors) 

F\a,P) = e-^ a P\(3,-a). (68) 

The quantization with respect to the finite frame {\a, P)} s a (9= _ s , we call finite frame 
quantization, is defined as the linear map / h-> Af, where 

f:{-s,-s+l,...,s}x{-s,-s+l,...,s} — ► C 

is a classical observable defined on the finite phase space [2T] 133] and 



A f = \ E f(<*,P)\a,P){a,P\ (69) 

a,/3=—s 

the corresponding linear operator [11]. For the linear operator corresponding to the 
function 



namely, 



q:{-s,-s + l,...,s}x{-s,-s + l,...,s}^R, q(a, P)=aJ% (70) 



(71) 



a,/3=— s 



the basis vectors are eigenvectors (see Figure [2]) 

/ — s 

A q \e n ) = \ n \e n ), where \ n = J^f £ «g 2 («+n). 
Therefore, A g admits the representation 

s 



(72) 



(73) 
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From the relation (1681) . it follows that the linear operator corresponding to 
p:{-s,-s+l,...,s}x{-s,-s+l,...,s} — >R, p(a,/3) = /3 1 f 
namely (see Figure EJ, 

A - I / 2vr 
d\ d 



1\/ 2 T E P\a,/3)(a,l3\ 

a,/3=~s 



satisfies the relations 



A r . 



F + A q F, 



13 



(74) 



(75) 



(76) 



where \e n ) = F + |e n ). The linear operators A q and A p can be regarded as alternative 
versions of the usual coordinate and momentum operators [151 EH [26], [28] 



Q = \/T E n\e n )(e r 

n=—s 



and 



P 



n\e n ){e n \ 



(77) 



The set (see Figure [3]) 

Ad = { (A n , A fc ) | n,ke{-s,-s+l,...,s-l,s}} (78) 
can be regarded as an alternative version of the finite phase space [6] 

V <i= {( n \ftt> k yftt) | n,ke{-s,-s+l,...,s-l,s}Y (79) 

They are usually described by using the set {— s, — s+1, s}x{— s, — s+1, s—1, s} as 
a mathematical model. Numerically, one can remark the tendency A n ~ n^j2n/d for d 
large enough. In Figure El the central values of A q are almost equal to the corresponding 
eigenvalues of Q. This central part of the spectrum of A q becomes larger and larger for 
increasing d. 



The linear operator corresponding to the function / = \(p 2 + q 2 ), that is, to 
/(a,/3) = *(a 2 + /3 2 ) 



30) 



Figure 3. Two alternative versions of the finite phase space (case d=21). 
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namely [5] 

= § E^=- s (« 2 +/3 2 )l«,/3)(«,/3| (81) 

is the Hamiltonian of a finite oscillator. Numerically, one can remark that the energy 
levels are non-degenerated (see Figure [2]) and have the tendency to become 1, 2, 3, ... 
for d — > oo. In Figure El the lower-lying eigenvalues of ii^2) are almost identical to 
the corresponding eigenvalues of A\ z \2. This part of the spectrum becomes larger and 
larger for increasing d. 



5. Definition of an exact DFrFT based on the finite frame quantization (Fq) 

By using the property 

F\a,(3) = e -^\(3,-a) (82) 
of the finite frame {\a, (3)} s a a = _ s we get the relation 



Any eigenfunction of Ai 



Ai 



2 -H? 2 ) J 



(83) 



(P 2 +9 2 ) 



is at the same time an eigenfunction of the finite Fourier 
transform F. The finite oscillator with the Hamiltonian Aiu^ 2 ) is a finite counterpart of 
the translated harmonic oscillator A\ z \2. The normalized eigenfunctions $ m of Ai^j^y 
considered in the increasing order of the number of sign alternations, represent a finite 
counterpart of the Hermite- Gaussian functions \l/ m . Since <3> m is an eigenfunction of F 
corresponding to the eigenvalue (— i) m we have 

' d-i 



.771=0 



ip{n). 



54) 



The transformation 



where 



H 



Q 



E 

m=0 



e --ma 



i/j(n) 



55) 



represents the definition we propose for the discrete fractional Fourier transform. The 
operator Fq is an a-th power of F. Our transformation Fq, as well as F^, is unitary 

(86) 



Qi 



I 



and satisfies the relations 

F°=I, F Q = F and Fq Fq — Fq +/3 . (87) 
The Harper functions h m correspond to the normalized eigenvectors of the matrix 



(2(cos ^-2)8 jk + 6 j)k+1 + 5 j)k _i + 5 jtk _2s + ^j,k+2 S )_ 



■s<j,k<s 
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Table 1. The values of the difference Ah (m) — Aq (m) in the case d = 21. 
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m 


A H (m)-A Q (m) 


m 


A H (m)-A Q (m) 


m 


A ff (m)-A Q (m) 





0.003234 


7 


0.085618 


14 


0.155817 


1 


0.008071 


8 


0.126788 


15 


0.067190 


2 


0.013610 


9 


0.102790 


16 


0.173158 


3 


0.022260 


10 


0.202624 


17 


0.036419 


4 


0.031979 


11 


0.114106 


18 


0.025582 


5 


0.047942 


12 


0.193413 


19 


-0.009772 


6 


0.060942 


13 


0.137640 


20 


0.002824 




Figure 4. The values of Ah (in) (bullets) and Aq(?ti) (squares) in the case d = 21. 



and our functions $ m correspond to the normalized eigenvectors of the matrix 



7T 

d 2 ^ 



cr 



e <i 



(a + k)g(a+j] 



a,f)=—s 



59) 



-s<j,k<s 

Our functions $ m approximate the Hermite-Gaussian functions much better 
than the Harper functions h m . By denoting 

A Q (m) = max |$ m (n) - ^/2-K/d^ m {nJ2n/d)\ 

~ s < n < s (90) 

A H (m) = max \h m (n) - {/2tt /d^ m (nJ2ir /d)\ 

—s<n<s 

we have (see Table [D and Figure HJ) 



Ag(m) < Ajy(m) for almost all mG{0, 1, 2, d— 1}. 
Therefore, the discrete fractional Fourier transform Fq[ 



(91) 



computed by using our 
definition (1851) approximates the corresponding continuous fractional Fourier transform 
better than the discrete fractional transformation F^[?/>] computed by using the 
definition ( )42|) based on the Harper functions. For example, in the case d = 21, the 
discrete counterpart of the Gaussian 



10 • 



W) 



-5x 2 



(92) 



New definition of the discrete fractional Fourier transform 
is the 21-periodic function (see Figure [5]) 

oo 



010 
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(93) 



and (see Figure [6]) 

n 5 [0io]W 
F^ 5 [Sio](n) 



but, 



and 



F°- 5 [f 10 }(n^/d), for ne{-s,-s+l, 
T°- 5 [f 10 \(ny^/d), for ne{-s,-s+l, 



8-1,8} 

s-l,s} 



max 

-s<n<s 

max 

-s<n<s 



max 

-s<n<s 

max 

-s<n<s 



max 

-s<n<s 

max 

-s<n<s 



^(F^ 5 [ gi0 ])(n) - mt(^[f 10 ])( nv ^/d)\ = 0.094066 



JKe(F°j 5 [g 10 ])(n) - mt(^[f l0 })(n^/d)\ = 0.0548102 



3m(F° H 5 [Q 10 ])(n) - 3m(^[f 10 })(ny^/d)\ = 0.0768065 
3m(F° Q 5 [ Ql0 ])(n) - 3m(F°- 5 [f 10 ])(n^/d)\ = 0.0496917 



(94) 



F^ 5 [flio]H -^°- 5 [/io](n>/27r/d)| = 0.0940663 



F°q 5 [0io]H - F°- 5 [fio](nV^/d)\ 



0.0567312. 



(95) 






Figure 5. Gaussian f w , its Fourier transform J"[/io] and the corresponding discrete 
counterparts gio and F[{jio]. 



6. Concluding remarks 

The finite frame {\a, fi)} s a ^ = _ s is a finite counterpart of the standard system of coherent 
states {|z)} 2e c. In terms of the Perelomov method, {|;z)} ze c can be obtained by using 
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Figure 6. The real part (left hand side) and the imaginary part (right hand side) of 
the 0.5-th order fractional Fourier transform of the Gaussian / 10 and of its discrete 
counterpart gio computed by using the definition based on the Harper functions (second 
line) and the proposed definition (third line). 

the Heisenberg-Weyl group [22]. The finite frame {\a, fi)} s a 8=-s can be generated in a 
very similar way by using a finite version of the Heisenberg-Weyl group [251 EHJ EE] . The 
Hamiltonian of the quantum harmonic oscillator can be obtained by using the coherent 
state quantization [TT]. The Hamiltonian of our finite oscillator is obtained in a very 
similar way by using the finite frame quantization [8] and its eigenfunctions approximate 
the Hermite-Gaussian functions better than the Harper functions. 

To our knowledge, the finite frame quantization has been presented for the first 
time in the papers of Gazeau et al. [71 HH 112] . Our finite oscillator, presented for 
the first time in [8], represents the main known application of this method. Discrete 
counterparts of the differential operators are usually obtained by using finite difference 
operators. Our results show that, in certain cases, the linear operators defined by using 
the finite frame quantization approximate the continuous operators better than those 
based on finite-differences. 

The finite frame quantization can be regarded as an extension of the method used 
to define Q and P. In all the cases the starting point is a resolution of the identity. By 
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following the analogy with the definitions of the coordinate and momentum operators 



I= E kn)( £ n| i-> Q=\Jlt E n\e n ){e n \ 

n=—s n=—s 

1= E \e n )(e n \ ^ P =\f^t E ™k~n>(£n| 

ri=— s v n=— s 



(96) 



we define the operator corresponding to a function / defined on the finite phase space 
1 = 2 E \<*,P)(a,l3\ ^ A f = \ E /(«,/?) (97) 

a,/3=—s a,f}=—s 

We have obtained an alternative definition of DFrFT more consistent with the 
continuous transform than the best known definition. At the same time, our results 
represent a significant contribution to the discrete version of the phase space quantum 
mechanics [271 EH [33] . 
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